home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dlae2.z / dlae2
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAAEEEE2222((((3333FFFF))))                                                            DDDDLLLLAAAAEEEE2222((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLAE2 - compute the eigenvalues of a 2-by-2 symmetric matrix  [ A B ]  [
  10.      B C ]
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE DLAE2( A, B, C, RT1, RT2 )
  14.  
  15.          DOUBLE        PRECISION A, B, C, RT1, RT2
  16.  
  17. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  18.      DLAE2  computes the eigenvalues of a 2-by-2 symmetric matrix
  19.         [  A   B  ]
  20.         [  B   C  ].  On return, RT1 is the eigenvalue of larger absolute
  21.      value, and RT2 is the eigenvalue of smaller absolute value.
  22.  
  23.  
  24. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  25.      A       (input) DOUBLE PRECISION
  26.              The (1,1) element of the 2-by-2 matrix.
  27.  
  28.      B       (input) DOUBLE PRECISION
  29.              The (1,2) and (2,1) elements of the 2-by-2 matrix.
  30.  
  31.      C       (input) DOUBLE PRECISION
  32.              The (2,2) element of the 2-by-2 matrix.
  33.  
  34.      RT1     (output) DOUBLE PRECISION
  35.              The eigenvalue of larger absolute value.
  36.  
  37.      RT2     (output) DOUBLE PRECISION
  38.              The eigenvalue of smaller absolute value.
  39.  
  40. FFFFUUUURRRRTTTTHHHHEEEERRRR DDDDEEEETTTTAAAAIIIILLLLSSSS
  41.      RT1 is accurate to a few ulps barring over/underflow.
  42.  
  43.      RT2 may be inaccurate if there is massive cancellation in the determinant
  44.      A*C-B*B; higher precision or correctly rounded or correctly truncated
  45.      arithmetic would be needed to compute RT2 accurately in all cases.
  46.  
  47.      Overflow is possible only if RT1 is within a factor of 5 of overflow.
  48.      Underflow is harmless if the input data is 0 or exceeds
  49.         underflow_threshold / macheps.
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.